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Abstract. The distribution patterns, sizes and nucleation sites of aggregated prisms on the outer 
shell surface were examined in 16 species of Bivalvia and modeled theoretically. Biometric analy- 
sis shows a negative correlation between the median size and variation of sizes of calcitic simple 
prisms. In species with aragonitic vertical composite prisms, instead, the density of prisms tends 
to decrease when their nucleation sites are randomly distributed. Comparison of the results of 
computer simulations with those of biometric analyses reveals the following: 1) a positive correla- 
tion between growth rate of prisms and the probability of nucleation for simple prisms, and 2) a 
limit of the number of nucleations per unit time in vertical composite prisms. Prism size correlates 
with the growth rate of the entire shell or prisms, and increases as the shell grows faster or prisms 
grow slower. 

Key words: biomineralization, bivalves, prismatic structure, shell growth rate, theoretical morpho- 
logy 



Introduction 

The microscopic features of a molluscan shell record cer- 
tain physiological conditions of the organism at the time the 
shell is formed. Carbonate minerals within a molluscan 
shell crystallize and grow under physicochemical conditions 
controlled by the physiology of the organism. Wada (1972, 
1985) reported a seasonal change of shape and size of 
aragonite crystals in the nacreous layer of such bivalves as 
Pinctada fucata, Pinna attenuata and Hyriopsis schlegeli. 
He suggested that the topography of growing crystals may 
reflect the rate of crystal growth or the degree of 
supersaturation of the extrapallial fluid. Quantitative analy- 
sis of the relationship between size and/or shape of crystals 
and the rate of crystal growth within a shell provides a reli- 
able basis to understand the ‘paleophysiology’ of fossil or- 
ganisms. 

The goal of this study is to clarify the relationship between 
the geometry of bivalve shell microstructure and the relative 
growth rate of crystals or of the entire shell. For this pur- 
pose, the present study focuses on the geometry of the 
outer shell surface of a simple or vertical composite pris- 
matic shell layer. Although they differ in the ultrastructure of 
prisms, those two prismatic structures both consist of many 
parallel-arrayed columnar units (Carter and Clark, 1985; 
Carter et al. , 1990). Each prism is surrounded and 
bounded by an organic matrix showing a honeycomb-like 
appearance on the outer shell surface. 



For understanding the rule or algorithm forming the geo- 
metric pattern of shell microstructure, theoretical morphol- 
ogy is particularly useful (Ubukata, 1997a, b, 2000). In the 
present study, a biometric analysis of size and nucleation 
sites of prisms in actual shells was carried out in 16 species. 
Furthermore, a theoretical morphological modeling of growth 
kinematics of aggregated prisms was attempted, and the 
computer simulations of that model were compared with the 
results of the biometric analyses. 

Biometric analyses 
Material and methods 

The outer shell surface of a simple prismatic or vertical 
composite prismatic outer shell layer was examined in 16 
extant species of Bivalvia (Table 1 ). Each species was rep- 
resented by a single specimen, except Anodonta woodiana. 
Most of them were collected at various localities around the 
Japanese Islands and the Philippines. All the specimens 
examined are stored at Shizuoka University (SUM). 

In order to remove the periostracum from the shell com- 
pletely, the shells examined were bleached for one day. 
Pieces of them were washed, dried in air, coated with gold 
using a JEOL JFC-1500 ion coater, and examined with a 
JEOL JSM 5800LV scanning electron microscope operated 
at 15kV and interfaced to a computer (Dell Optiplex Gxa 
EM). 

In order to analyze size-frequency distribution of prisms, 
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Table 1. List of material examined. All specimens have the prefix SUM. 



Family 


species 


locality 


specimens 


Pteriidae 


Pteria penguin (Roding) 


San Luice, Bathangas, Philippines 


HM-B-0014 




Pinctada maculata (Gould) 


Iriomote Is., Okinawa, southwest Japan 


HM-B-0015 


Isognomonidae 


Isognomon perna (Linnaeus) 


Iriomote Is., Okinawa, southwest Japan 


HM-B-0016 




/. ephippium (Linnaeus) 


San Luice, Bathangas, Philippines 


HM-B-0017 


Malleidae 


Malleus regula (Forsk&l) 


Iriomote Is., Okinawa, southwest Japan 


HM-B-0018 


Pinnidae 


Atrina pectinata (Linnaeus) 


Ariake, Saga, western Japan 


HM-B-0019 




A. vexillum (Born) 


Honda Bay, Palawan, Philippines 


HM-B-0020 


Ostreidae 


Crassostrea gigas (Thun berg) 


Misaki, Kanagawa, Central Japan 


HM-B-0021 


Margaritiferidae 


Margaritifera laevis (Haas) 


Nakagawa, Hokkaido, northern Japan 


HM-B-0022 


Unionidae 


Inversidens reiniana (Kobelt) 


Lake Biwa, Shiga, Central Japan 


HM-B-0023 




Unio biwae Kobelt 


Lake Biwa, Shiga, Central Japan 


HM-B-0024 




Lanceolaria oxyrhyncha (Martens) 


Lake Biwa, Shiga, Central Japan 


HM-B-0025 




Anodonta woodiana (Lea) 


Lake Biwa, Shiga, Central Japan 


HM-B-0026, -0027 




A. calypygos Kobelt 


Lake Biwa, Shiga, Central Japan 


HM-B-0028 




Cristaria plicata (Leach) 


Lake Biwa, Shiga, Central Japan 


HM-B-0029 


Trigoniidae 


Neotrigonia margaritacea (Lamarck) 


French Is., Australia 


HM-B-0030 





areas of prisms on the outer shell surface were measured at 
6-12 positions along a growth increment on the shell sur- 
face. An SEM image of the measured portion was saved as 
a computer bitmap file (Figure 1A). Next, the boundaries 
between prisms were traced on a NEC PC-9821 VI 66 per- 
sonal computer using Microsoft PowerPoint 7.0, and then 
each prism was colored differently using Justsystem Hanako 
PhotoRetouch (Figure IB). Subsequently, the area of each 
prism was measured by counting pixels. For this counting, 
a program written in VISUAL BASIC was used on a personal 
computer. 

Since the size-frequency distribution of the areas of 
prisms is generally right-skewed, the mean and standard de- 
viation are not suitable for representing the distribution of the 
areas of prisms. Therefore, the “average” area of prisms on 
a shell is represented by the median of the areas (5), and the 
variation of the areas is expressed as a standardized hinge 
spread (Q x ), which is defined as follows: 



where Qa nd Q, are the first and third quartiles of the areas 
of prisms, respectively (Hoel, 1976). S and Q s were both 
estimated in all shell portions examined. 

On the outer surface of a vertical composite prismatic 
shell layer, microgrowth increments are clearly visible within 
a bleached prism (Figure 2A). In a simple prismatic shell 
layer, growth increments within a prism are faintly observed 
on the bleached shell surface (Figure 2B). In either case, 



Figure 1. A. SEM photograph of the outer shell surface of the 
vertical composite prismatic layer in Anodonta woodiana (SUM- 
HM-B0027), scale: 50 pm. B. Trace of the outlines of prisms on 
the SEM image of A. 
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Figure 2. A. SEM photograph of the outer shell surface of the vertical composite prismatic shell layer in A. woodiana (SUM-HM-B0027) 
showing clear growth increments within individual prisms, scale: 20 pm. B. The simple prismatic layer in Pinctada maculata (SUM-HM- 
B0015). Arrows indicate nucleation sites, scale: 20 pm. C. The SEM image of A. woodiana subdivided into squares, scale: 50 pm. D. 
Distribution of nucleation sites of prisms in C. 



the center of the circular growth increments is regarded as 
the nucleation site of the prism. 

The distribution of nucleation sites of prisms was analyzed 
in all shell portions in which the areas of prisms were meas- 
ured. An SEM image was divided into a squared grid of ap- 
propriate size for each square to include an average of two 
nucleation sites (Figure 2C). Next, the number of nuclea- 
tion sites was counted in each square (Figure 2D). The 
numbers of squares and nucleation sites in each square 
both define the Is index of Morisita (1959), which represents 
the nonuniformity of a distribution independent of the size of 
the quadrates. The Is index is defined as: 



1 * 



\x,(x, -1) 



( 2 ) 



where n is the number of squares, x, is the number of sites 
in the square /. The value of I& is zero wh$n the distribution 
is perfectly uniform (Figure 3A). h increases as the distribu- 
tion becomes nonuniform (Figure 3B), and it is expected to 
be one when sites are distributed randomly (Figure 3C). 
When the distribution of sites is biased considerably, Is has 
a large value (Figure 3D). The value of Is was estimated in 
every shell position examined, and the nonuniformity of the 
distribution of nucleation sites was represented by Is. 

Results 

The biometric analyses revealed a negative correlation 
between S and Q, in a single specimen of Isognomon 
perna, and a positive correlation in a single specimen of 
Malleus regula (Figure 4B). In another species, no correla- 
tion was observed, though the number of measured portions 
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Figure 3. Effect of the spatial distribution pattern on the value of Is. In this figure, the area is divided into 16 contiguous squares each 
including a number of or no dots. A. Perfectly uniform distribution. B. Nearly uniform distribution. C. Random distribution. D. 
Considerably concentrated distribution. 



within a single valve was not large enough to be significant. 
On the other hand, when the data from all species are com- 
bined, Qs is negatively correlated with S at the 0.01 level of 
significance, both for simple prismatic and vertical composite 
prismatic shells (Figure 4A-B). This shows that variation in 
the size of prisms tends to decrease as prism size in- 
creases. The negative correlation between Q s and S is 
clear especially in species with a simple prisms (Figure 4B), 
while the correlation is more or less obscure in species with 
vertical composite prisms (Figure 4A). 

A positive correlation between h and Q s is found only in 
each specimen of Unio biwae, Anodonta woodiana and 
Malleus regula (Figure 4C, D). However, when all the data 
are combined, a positive correlation clearly emerges be- 
tween h and Q, in species with vertical composite prisms 
(Figure 4C). In the case of species with simple prisms, Q s 
is positively correlated with h at the 0.05 level of signifi- 
cance, though the correlation is graphically unclear (Figure 
4D). The positive relationship between h and 0, indicates 
that variation of the size of prisms tends to increase with in- 
creasing randomness of the distribution of nucleation sites. 

Theoretical morphology 
Growing circles model 

To better understand the relationships between geometric 
features of prisms and the growth rate of each prism and/or 
of the net growth rate of the entire shell, the growth of aggre- 
gated prisms was modeled theoretically. 



For modeling the process of microscopic growth, observ- 
ing the initial growth stage of prisms helps understand the 
nature of crystal growth. Formerly, I reported that many 
small hemispherical incipient prisms occur on the inner sur- 
face of the periostracum at the growing margin in species 
possessing simple and vertical composite prisms (Ubukata, 
1994, pi. 2, figs. 1-3). Consequently, a growing-circles 
model, which represents the growth of aggregated prisms, is 
introduced here. 

Growth of a prismatic shell layer consists of three ele- 
ments, namely, nucleation of prisms, growth of prisms, and 
accretionary growth of the entire shell. During a single 
short growth step, the mantle secretes calcium carbonate 
and nucleation of prisms occurs within the nucleation zone 
on the inner surface of a periostracum (Figure 5A). The 
periostracum subsequently secreted by the mantle edge 
pushes the earlier produced periostracum and its incipient 
prisms into a more proximal part of the shell (Saleuddin and 
Petit, 1983). After the prisms pass through the nucleation 
zone, prisms gradually grow and elongate, forming a colum- 
nar structure. 

Let us consider hypothetical shell growth (Figure 5B). 
The accretionary growth of the entire shell during a short pe- 
riod of time is reflected in a shift of the nucleation zone. 
Growth of a prism is represented by the kinematics of an en- 
larging circle. Potential nucleation sites are distributed uni- 
formly within a nucleation zone of width h. The distance 
between the potential nucleation sites can be expressed by 
d, which represents the size of a unit cell. Growth compo- 
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Figure 4. A-B. Relationship between S and Q , in actual shells. A. Data from species with vertical composite prisms. B. Data from 
species with simple prisms. C D. Relationship between U and Q . for species with vertical composite prisms (C) and simple prisms (D). 



nents are expressed as functions of growth stage s, rather 
than as functions of time, since the time scale of the growth 
process is difficult to ascertain in many cases. Over a pe- 
riod of one growth step, the mantle secretes a periostracum 
at the shell margin, giving rise to the stippled area in Figure 
5B. If the growing margin of the shell shifts downward by d 



during the step, the growth step is regarded to be a unit in- 
terval of the growth. Then, the growth step As is generally 
defined as: 

AS = ~. ( 3 ) 

a 
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Figure 5. A. Schematic diagram of the radial section of the shell and mantle margin of a bivalve showing the position where nucleation 
of prisms occurs. B. The growing circles model. Black bold points indicate initiation sites of produced prisms, and gray ones potential nu- 
cleation sites of unborn prisms. During a given growth step As, the shell margin shifts by Al (stippled area) and the radius of prisms increases 
by A R (shaded area). The dotted line indicates the dorsal limit of the nucleation zone. 



where the growing margin of the shell shifts by Al during the 
growth step. 

Meanwhile, nucleation of a prism occurs at a potential nu- 
cleation site within the nucleation zone with a probability of 
q. Each prism is approximated by a circle which enlarges at 
a steady rate. During a growth step As, calcium carbonate 
precipitates along the circumference of each prism giving 
rise to a new additional rim shown by the shaded portion in 
Figure 5B, and the radius of each circle increases by A R. 
As the prisms grow, neighboring prisms come closer and fi- 
nally in contact with one another, as a result forming a 
boundary between two prisms. Nucleation of prisms occurs 
randomly during each growth step, as a result of irregularity 
of the settling time among prisms. Consequently, a growing 
circle often occupies the space of nucleation and/or growth 
of a neighboring ‘unborn’ prism. A newborn prism some- 
times loses in competition for space between neighboring 
prisms, and is geometrically terminated (Grigor’ev, 1965). 

In the growing-circles model, growth of a shell and prisms 
is generally expressed by the following three parameters: C 



: the standardized growth rate of prisms, defined as the in- 
crease of the radius of a prism per growth step, normalized 
by size of a unit cell (d), P : probability of nucleation per 
growth step in each potential site, and L: the extent of the 
nucleation zone standardized by d. Probability of ‘failure’ of 
nucleation at each site per growth step is expressed as 1- 
P, and the probability of failure of nucleation during a growth 
step As (=1-<?) is obtained by raising 1 -P to As th power. 
Then, three growth parameters C, P and L are given by: 

Lmh-. (4) 

Al d 

Now, we can generally define a growth increment during 
an arbitrary growth step if three parameters C, P and L are 
given. 

Computer simulation 

In order to evaluate the effects of parameters C, P and L 
on the geometric pattern of prismatic structure, computer 
simulations were performed for growth of prisms based on 



increase /, 
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Figure 6 Color display of the growing circles model. Each prism is identified by its color. As C increases, size variation of the prisms 
also increases As P or L increases, the median size of prisms tends to decrease. 
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Figure 7. Three-dimensional block diagrams showing the relationships among growth parameters (C, P, L) and statistics S (A-C), Q< 
(D-F) or h (G-l) on a C-P diagram at L values of 1 .2, 3.6 and 7.2. 



the growing-circles model. In the theoretical model, the 
area of a hypothetical prism was measured as the number of 
pixels (A) on the display surface, and both the standardized 
hinge spread ( Q x ) and the median of a standardized area of 
hypothetical prisms (S sl ) defined beiow were estimated in 
each model: 

s.,=4- ® 

a 



where A is the median of A. The dimension of d is ex- 
pressed as pixels on the computer. The coordinates of a 
nucleation site were recorded on each hypothetical prism for 
calculating the value of h in a model. Computer simulations 
were carried out with a program written in VISUAL BASIC by 
means of a 64-bit workstation computer (Visual Technology 
VT-Alpha 600) interfaced with a CRT (liyama A702H). 

Figure 6 shows a spectrum of geometric patterns of hypo- 
thetical prisms that were made by the growing circles model. 
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A 



B 




Figure 8. Summary of the biometric analyses and computer simulation. From coupling their results with the consequences from defi- 
nition of the parameters, a relationship between growth parameters is inferred either in simple prisms (A) or composite prisms (B). 



Each prism is identified by its color in a theoretical prismatic 
structure. If C is large, each prism grows considerably while 
it passes through the nucleation zone, as a result of the 
great irregularity of the birth time among prisms. Then, as 
C increases, the variation of prism size also increases. As 
P or L increases, the number of prisms increases while their 
size decreases. 

Three-dimensional scatter diagrams illustrating the values 
of Su in relation to C and P, at L= 1.2, 3.6 and 7.2 are given 
in Figure 7A-C, respectively. which means an ‘average’ 
prism size, seems to be inversely proportional to P, which 
represents the probability of nucleation. In addition, at the 
same combination of C and P, S., tends to increase as L de- 
creases. This fact indicates that the size of hypothetical 
prisms increases as the nucleation zone becomes narrow. 

Figure 7D-F is a series of diagrams showing the relation- 
ship among four parameters C, P, Q and L. When L is large 
enough, Q, is positively correlated with C (Figure 7E, F). 
This fact suggests that variation of prism size tends to in- 
crease with increasing growth rate of prisms (A R in Eq.1) or 
with decreasing the accretion rate of the entire shell (A/ in 
Eq.1). When L is small, no clear relationship exists be- 
tween pairs of C, P and Q, (Figure 7D). In a high C condi- 
tion, Q, at a given P also tends to increase with increasing L 
(Figure 7D-F). 

Figure 7G-I shows the variation of h in relation to C, P and 
L. When L is large, h gradually decreases with increasing 
P. Since h represents the degree of nonuniformity, this re- 
sult indicates that the nucleation sites tend to be distributed 
uniformly as the probability of nucleation increases. In the 
high P region, h at a given C decreases as L increases. 
This fact indicates that nucleation sites tend to be distributed 
randomly as the nucleation zone becomes narrower, when 
the probability of nucleation is high enough. 



Geometric pattern and growth rate of prismatic shell 

As mentioned above, the biometric analyses indicated a 
negative correlation between S and Q< especially in the 
species with simple prisms (Figure 4B), and the computer 
simulation predicted an inverse relationship between P and 
S„ (Figure 7A-C) and also a positive correlation between C 
and Q s (Figure 7E, F). To sum up these results, it can be 
predicted that P increases as C increases for those species 
with simple prisms (Figure 8A). That is also inferred from 
the definition of parameters. Since both P and C are in- 
versely proportional to the growth rate of the entire shell, a 
positive correlation between C and P is quite reasonable 
(Figure 8A) if the size-frequency distribution of prisms is con- 
trolled mainly by the growth rate of the entire shell. 
Furthermore, even if both the growth rate of prisms and the 
probability of nucleation reflect the activity of mantle secre- 
tion, and if the size of prisms is controlled mainly by the ac- 
tivity of mantle secretion, the positive correlation between C 
and P is also expected. 

Coupling of the biometrics and the simulation also sug- 
gests to us a relationship between the parameters P and L. 
The biometric analyses demonstrated a positive correlation 
between Is and £?,, especially for species with vertical com- 
posite prisms as shown in Figure 4C. In addition, the com- 
puter simulation predicted a positive relationship between L 
and Q., and a negative correlation between P and h as 
shown in Figure 7G-I. These two results suggest that P 
possibly decreases as L increases in species with vertical 
composite prisms (Figure 8B). Actually, an inverse relation- 
ship between L and P is expected if we assume that the fre- 
quency of nucleation per unit time interval is fixed. 

To ascertain this prediction, relationships among simu- 
lated values of S,„ h and Q, were analyzed for each combina- 
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a a 







Figure 9. Relationships between 5„ and & and between / 5 and g, when three parameters satisfy equations as follows: C-aP , P=b!L , 
where a and 6 are coefficients. N represents the total number of samples examined, r exhibits the correlation coefficient, and p is the signifi- 
cance level of the correlation. A-B. a- 3, 0.05 < b <; 0.175. C-D. 1 <; a <; 3, 6=0.05. 



tion of C, P and L, when a positive correlation between C and 
P and an inverse relationship between L and P are assumed, 
as in the following relations: 

C = aP, P = (6) 

where a and b are coefficients. Computer simulations were 
performed under the following conditions: 1) a is fixed to 3, 
and b varies between 0.05 and 0.175, 2) a varies between 1 
and 3, and b is fixed to 0.05. Figure 9 is a diagram showing 
the relationships between s M and Q s and between / 5 and Q s . 

The results of the simulations under a fixed value of a (=3) 
and various values of b show a negative correlation between 



S sl and Q s (p<0.01, Figure 9A), but do not indicate any linear 
relationship between h and Q s (Figure 9B). Under even this 
condition, a negative correlation between S sl and Q s and a 
positive correlation between I 6 and Q s were approached by a 
multiple regression analysis, which shows a significant trend 
for Q s to increase as a function of - S s , and h: 

Q; =-0.5735/ + 0.128 V (r= 0.565, F=11.94, p<0.01), 

where Q s * , S sl * and h* are standardized variables of Q s , 
S st and Is, respectively. In this case, a negative relationship 
between S s , and Q s is especially prominent. This relation- 
ship, concordant with that of the biometric analyses as 
shown in Figure 4B, D, strongly supports conspicuous posi- 
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tive correlation between C and P assumed as in the Eq. 6, 
particularly for simple prisms. As stated above, a positive 
correlation between C and P suggests that both prism size 
and its variation are mainly controlled either by the growth 
rate of the entire shell, or by the activity of mantle secretion. 
If the former is the case, simple prisms must tend to be uni- 
formly large as the entire shell grows faster. On the other 
hand, if the latter is the case, density of nucleation and 
growth rate of prisms both must decrease as the secretive 
activity of the mantle decreases, and a negative relationship 
is expected between size and the growth rate of simple 
prisms. 

On the other hand, the results of computer simulations 
with various values of the coefficient a and the fixed coeffi- 
cient b (=0.05) show a positive correlation between h and 
Q< (p<0.01, Figure 9D), but no significant relationship be- 
tween Sit and Q, (Figure 9C). A multiple regression analysis 
for Q ( on S v and Is provides a significant trend to increasing 
Q as a function of S r and h: 

Q s m =-0.3935/ +0.686 Is' {>- 0.582, F=9.99, /?<0.01). 

In this equation, a positive relationship between Is and Q> 
is more striking than a negative relationship between and 
Q t . This is concordant with the results of the biometric 
analyses in species with vertical composite prisms shown in 
Figure 4A, C. This fact supports the assumption that P has 
an inverse relationship with L as defined in Eq. 6, particularly 
in vertical composite prisms. This assumption implies that 
the number of nucleations per unit time interval is fixed to a 
constant value. In the species having a large nucleation 
zone, in which the value of L is large, maintenance of the 
probability of nucleation inevitably causes a large number of 
nuclei. If we assume an upper limit of the total number of 
nuclei per unit time interval, the probability of nucleation will 
vary inversely with the width of the nucleation zone. If the 
total number of nuclei per unit time interval is fixed, the den- 
sity of nucleation is expected to decrease as the total shell 
grows faster. In this case, the median size of prisms is con- 
trolled mainly by the growth rate of the entire shell. This 
fact suggests that, in the case of vertical composite prisms, 
prism size tends to increase as the entire shell grows faster. 

Discussion and conclusion 

Wada (1961) studied the size of crystals in a nacreous 
shell layer of Pinctada martensii in relation to the rate of cal- 
cium carbonate deposition. He demonstrated that when the 
rate of shell deposition is at a maximum, a large number of 
small crystals occurs on all the nacreous surfaces, while 
larger crystals occur as the rate of deposition decreases. In 
addition, Wada (1972, 1985) also reported an inverse 
relationship between them for such bivalves as Pinctada 
fucata , Pinna attenuata and Hyriopsis schlegeii . Wada as- 
sumed that the primary factor determining size of crystals 
was the degree of calcium carbonate concentration in the 
extrapallial fluid. He regarded the rate of calcium carbonate 
deposition as the rate of crystal growth, and thought that 
larger crystals tend to grow slowly at a low degree of 
supersaturation of the fluid, at which the frequency of nu- 
cleation diminishes. The results of the present study may 



partly support those of Wada, since the negative correlation 
between size and growth rate of crystals is also expected in 
this study if the size of prisms is assumed to be controlled by 
the activity of mantle secretion. 

Unlike prismatic structure, nacreous structure does not 
form the outermost shell layer in bivalves. In the nacreous 
layer, deposition corresponds to thickening of the shell or 
growth of crystals, rather than growth of the entire shell. 
Therefore, the size of crystals in the nacreous layer corre- 
lates with the rate of crystal growth rather than with the 
growth rate of the entire shell. On the contrary, the net 
growth rate of the entire shell, which reflects the growth rate 
of the soft parts, seems to be significant for the size of crys- 
tals which constitute the outermost shell layer. 

Ubukata (1994) claimed that relatively rapid growth of 
prisms produces prisms prominently inclined to the outer 
shell surface, because of retardation of the initiation of their 
forward growth relative to the radial direction. Although 
such inclined prisms are commonly found in species belong- 
ing to Unionidae and Ostreidae (Ubukata, 1994), their 
prisms also characteristically fell in low (small size) and high 
Q< (irregular in size) regions in Figure 4A, B (refer to Table 1 
for higher taxonomy). The present study suggests that 
such a pattern is produced under the condition of either 
rapid growth of prisms or slow growth of the entire shell. 
Ubukata (1994) demonstrated that the relative growth rate of 
prisms to that of the entire shell determines the orientation of 
elongation of the prisms. 

This study suggests that the sizes of prisms and their size 
variability correlate with the growth rate of the shell in bi- 
valves. It is well known that internal microgrowth incre- 
ments often provide a high-resolution record of growth rate 
and physiological condition of a bivalve (Lutz and Rhoads, 
1980), particularly in species with crossed lamellar and hori- 
zontal composite prismatic shells. In simple prismatic and 
vertical composite prismatic shells, however, it is hard to es- 
timate the growth rate of the shell from an analysis of inter- 
nal microgrowth increments, because periodic growth incre- 
ments tend to be obscure. The size-frequency distribution 
of prisms may be a promising index for estimating the 
growth rate of a fossil bivalve shell. 

This study indirectly estimates the relationship between 
the size of prisms and the growth rate of the shell or prisms 
from computer simulations, but the direct observation of the 
growth rate of the shell or prisms remains to be performed. 
An experimental study on growth of a shell and/or crystals 
may be required for clarifying the relationship between the 
crystal size and the growth rate more clearly, and for obtain- 
ing a regression equation to predict the growth rate from 
prism size. The growth rate of the shell and/or crystals is 
clearly related to the physiological condition of a fossil or- 
ganism, and has the potential to in spire a research field of 
‘paleophysiology’. 
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